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Bose-Einstein Condensation of Magnons in TlCuCl 3 : Phase diagram and specific heat 
from a self-consistent Hartee-Fock calculation with a realistic dispersion relation. 

Gregoire MiSGUiCH 1 and Masaki Oshikawa 2 

Service de Physique Theorique, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France 
2 Department of Physics, Tokyo Institute of Technology, Oh-oka-yama, Meguro-ku, Tokyo 152-8551 Japan 

We extend the self-consistent Hartree-Fock-Popov calculations by Nikuni et al. [Phys. Rev. 
Lett. 84, 5868 (2000)] concerning the Bose-Eistein condensation of magnons in TICUCI3 to 
include a realistic dispersion of the excitations. The result for the critical field as a function 
of temperature behaves as H C (T) — H c (0) ~ T 3//2 below 2 K but deviates from this simple 
power-law at higher temperature and is in very good agreement with the experimental results. 
The specific heat is computed as a function of temperature for different values of the magnetic 
field. It shows a A-like shape at the transition and is in good qualitative agreement with the 
results of Oosawa et al. [Phys. Rev. B 63, 134416 (2001)]. 



1. Introduction 

TIC11CI3 is a spin-i magnetic insulator with a spin 
gap 1 of A = 7.5 K. 2 It has been successfully described 
as copper dimcrs with an intra-dimer antiferromagnetic 
exchange energy J ~ 5.5meV and weaker (< 1.5meV) 
inter dimer couplings. 3,4 At zero temperature, an ap- 
plied magnetic field H closes the gap at the critical 
field H = H c (0), giving rise to a quantum phase tran- 
sition. H c (0) is related to the gap by g^BH c (0) = A, 
where [Ib is the Bohr magneton and g is the Lande g- 
factor. The field-induced phase transition continues to 
finite temperature T, with the temperature-dependent 
critical field H C (T). Above the critical field, a magnetic 
long-ranged order in the plane perpendicular to applied 
field develops. 5 ' 6 The existence of the ordering transition 
was predicted by a standard mean-field theory for spins. 7 
However, several characteristic features of the transition 
could not be explained by the mean-field theory. The 
two most notable features are the cusp-like minimum of 
the magnetization as a function of the temperature at 
the transition, and the power-law like dependence of the 
critical field 

H C (T) - H c (0) cx X"* (1) 

in the low temperature regime. The mean-field theory 7 
rather predicts a monotonic decrease of the magnetiza- 
tion and an exponentially fast approach of the critical 
field H C (T) to its zero-temperature limit H c (0), on low- 
ering the temperature. 

These features were successfully explained, at least 
qualitatively, as a Bose-Einstein condensation (BEC) of 
spin triplet excitations (magnons). 8 The cusp- like mini- 
mum of the magnetization at the transition temperature 
is understood with the decrease of the non-condensed 
magnons at all temperatures and the increase of the con- 
densed magnons below the transition, as the temperature 
is lowered. Moreover, the self-consistent Hartree-Fock- 
Popov (HFP) approximation on the magnon condensa- 
tion gives the power-law dependence (1) with the expo- 
nent <f) = 3/2, if the dispersion of the magnons is taken 



to be quadratic. 

As one can easily control the magnetic field, which 
corresponds to the chemical potential of the magnons, 
this system provides a new arena for the study of BEC, 
in a grand-canonical ensemble with a tunable chemical 
potential. 9 

However, the results of the HFP approximation given 
in Ref. 8 are not quite satisfactory to describe the experi- 
mental data in a quantitative manner. In order to further 
extend the study of magnon BEC, it would be important 
to improve the HFP approximation and clarify its range 
of validity. 

One of the problems is that the HFP approximation 
predicts a discontinuous jump of the magnetization at 
the transition temperature, which is not observed. This 
is considered to be an artifact of the HFP approximation, 
and related to its breakdown due to strong fluctuation 
in the vicinity of the transition. In this paper we rather 
focus on another problem concerning the phase bound- 
ary. That is, while the experimental results are roughly 
in agreement with the power law (1), the reported val- 
ues 5, 6,8,10,11 of the exponent <j> = 1.67 ~ 2.2 are con- 
sistently larger than the HFP prediction 3/2. Although 
it was suggested that the deviation is again due to the 
fluctuation effects, it has not been clarified. 

In the present work, we extend the self-consistent HFP 
calculations 8 by including a realistic dispersion calcu- 
lated from microscopic models 3 ' 4 instead of the quadratic 
approximation ~ w— used previously. 8 The critical 
field H C (T) obtained by this method is in very good 
agreement with the experiments and represents a signifi- 
cant improvement over the simple quadratic approxima- 
tion. Therefore the puzzle regarding the discrepancy of 
the exponent <f> between the theory and the experiment 
is solved within the HFP framework. Here we note that 
there are related theoretical works 12-15 on this problem. 
We will comment on them later in Discussions. 

We also make several other checks of the HFP approx- 
imation with the experimental data, to show that HFP 
framework has a rather wide range of validity but the 
quadratic approximation fails above a rather low tern- 
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peraturc <~ 1 K for TICUCI3. Finally, the specific heat is 
also computed and compared with the results of Oosawa 
et al. 10 

2. Hamiltonian 

As in Ref. 8, the Zeeman splitting is assumed to be suf- 
ficiently large compared to temperature so that only the 
singlet and the lowest triplet states of each dimer need 
to be considered. With this approximation the system is 
described by an hard-core boson Hamiltonian 



Ti. — Tin + T~tu 



(2) 



TLk contains the zero-temperature magnon dispersion re- 
lation ek + A and the external magnetic field H: 

Kk = $>£&fc(e fc -/i) (3) 
k 

H = 9HbH — A (4) 

where it is assumed that eo = 0. The magnon-magnon 
interactions are described by 



Hu = ^N U i b i b l b k+i b k'- 

q.k,k' 



(5) 



and we will neglect the g-dependcncc of U q and set U q = 
U. As discussed by Nikuni et al., 8 this system undergoes 
a phase transition between a normal phase (at low field 
or high temperature) where the system is populated by 
thermally excited triplets to a "supcrfluid" phase where 
the bosons condense. This condensation is equivalent, in 
the spin language, to a field-induced three-dimensional 
magnetic ordering. 

3. Hartree-Fock-Popov treatment of the con- 
densed phase 

We reproduce the Hartree-Fock-Popov (HFP) mean- 
field analysis of Eqs. 2-5 which was discussed in Ref. 8. 
For a strong enough magnetic field the zero-momentum 
state (we assume that e k has a single minimum at k = 
0) is macroscopically occupied b\. =a = b k= o = \/N c = 
y/Nn c . From this we can write the interaction part of the 
Hamiltonian in terms of a constant, 2-, 3- and 4-boson 
operators: 16 



TLu — 

Ho = 

H 2 = 

n 3 = 

H4 = 



Hq + H 2 + H 3 + H4 
1 



— U iV r 2 
2N 



UNc 
N 



U- 



1 

h . q 



N 



b l) 



H.c 



2b\ b q 



U 
2N 



>b l b l b k + q bk' 

q.k,k' 



(6) 
(7) 

(8) 
(9) 
(10) 



where J2' means that the terms with creation or anni- 
hilation operators at k — are excluded. We perform 
a simple mean-field decoupling for TLu- While H3 gives 



zero in this approximation, H4 gives: 
Hf F = -U N{n - n c f + 2(n - n c )U ^ ' b { b k (11) 

k 

where n is the total boson density; it must be determined 
self-consistently from the thermal average over the spec- 
trum of n MF = n a + n 2 + nf F : 



n 



MF 



^+E' e ~ fe b l bk 



+ ^E'(^^ + & -, & ( 12 ) 



C 



,eff 



.off 



fj, 



fj, - 2Un 



UN 



-n\ -(n- n c f 



(13) 
(14) 



The mean-field Hamiltonian in the normal phase is ob- 
tained by setting n c = in the previous expression (al- 
ready in a diagonal form). In that case, the self-consistent 
equation for the density is 



n 



d 3 k 



/s(efc) 



(15) 



where Jb{E) = l/(exp(/3_E) — 1) is the Bose occupation 
number. 

When n c > 0, H MF can be diagonalized by the stan- 
dard Bogoliubov transformation: 



n 



MF 

Ek = 

b k = 

u k = 



J2'E k (a{ a k + 1 -) 

k 

\Je k 2 - (U n c f 
u k a k - v k aL fc 



i^'e~ fe + C(16) 



h 1 
2E k 2 



v k = 



2E k 



(17) 
(18) 

(19) 



The existence of a condensate (n c > 0) is possible when 
Ek is gapless, which implies /x cff = —Un c , or equiva- 
lent ly: 

g/j, B H = A + U (2n - n c ) (20) 

and n c are thus linearly related in the condensed 
phase and the self-consistent equation is now: 8 



n — n c = 



d 3 k 
J2wf 



^(f B (E k ) + l) 



(21) 



4. Dispersion relation for TICUCI3 

The dispersion relation of triplet excitations in 
TlCuCl 3 was measured at T — 1.5 K with inelastic neu- 
trons scattering by Cavadini et al. 3 This dispersion re- 
lation was very well reproduced by Matsumoto et al. 
within a bond-operator formalism. Their result is: 

e fe _ feo +A = y/(J + a k )*-al (22) 
a k = J a cos(k x ) + J aic cos(2k x + k z ) 

+2J abc cos(k x + k z /2) cos(fc y /2) (23) 
J = 5.501 meV , J a = -0.215 mcV(24) 
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Tabic I. Estimations of the gap (or critical field at zero temper- 
ature) from experiments. Mag. stands for magnetization, INS 
for inelastic neutron scattering, ESR for electron spin resonance, 
ENS for elastic neutron scattering (observation of the magnetic 
ordering) and C v for specific heat measurements. 



Rcf. 


A 


Method 


Shiramura et al. 2 


7.5 K 


Mag. 


(1997) 


(jHc = 5.6T) 




Tanaka et a/. 18 


7.68 K 


ESR 


(1998) 


(160 GHz) 




Oosawa et al. 


7.54 K 


Mag. 


(1999) 


(f H c = 5.61T) 




Tanaka et al.° 


7.66 K 


ENS 


(2001) 


(§if c = 5.7T) 




Cavadini et oi. a 


9.28K 


INS 


(2001) 


(0.8meV) 




Oosawa et al w 


7.66K 


C v 


(2001) 


(§ff c = 5.7T) 




Oosawa et al. u> 


7.54K 


INS 


(2002) 


(0.65mcV) 




Riicgg et al. 2u 


8.2K 


INS 


(2003) 


(0.71meV) 




Shindo et al LL 


7.33K 


Cv 


(2003) 


(f H c = 5.46T) 





Ja2c = -1.581 meV, J abc = 0.455 meV(25) 

where the Brillouin zone is doubled in the z direction 
(— 2ir < k z < 2ir) to represent the two magnon branches. 
The momentum shift by fc = (0, 0, 2ir) just insures the 
consistency between our convention that e — and 
the location of the minimum of the dispersion at fco in 
Refs. 3, 4 The dispersion relation above has a gap of 
0.7 meV, which is in agreement with the result of Ref. 3. 
However the studies based on a determination critical 
field as a function of temperature (see Table I) provide 
slightly smaller estimates for the gap (A ~ 0.65 meV) 
in TICUCI3. Therefore we corrected the value of J so that 
the dispersion relation is consistent with these data. The 
corrected value was chosen to insure A = 0.65 meV (or 
cquivalcntly (g/2)H c (Q) = 5.61T) : 

J = 5.489 meV (26) 

From the computation of curvature of around k = 
the effective inverse mass 17 1/m is 43.66 K (in units 
where ti 2 jks = 1), in agreement with the value taken 
in Rcf. 8. Fig. 1 shows the experimental data of Cava- 
dini et al. with the ej, given by Eqs. 22-26. The dotted 
line corresponds to the quadratic approximation; it only 
matches the full expression at very low energy. 

5. Critical density 

Within the HFP approximation the boson density n CT 
(or magnetization) at the transition is independent of 
the strength U of the magnon-magnon interaction as well 
as independent of the value of the zero- field gap Ao. It 
is obtained by setting /i cff = in Eq. 15. 21 If the full 
dispersion relation is used, the result has no adjustable 
parameter left. The result is shown Fig. 2 and is in good 
agreement with the experimental data. We note however 
that the discrepancy is larger when the field is applied 
along the b direction. We do not know the reason of the 




(0,27T,27r) (0,0, 2tt) (0,0,0) (tt,0,0) (0,0, 2tt) 



Fig. 1. Dispersion relation of triplet excitations. Full lines: Result 
of Eq. 22 with J given by Eq. 26. Dotted line: (anisotropic) 
quadratic approximation in the vicinity of the minimum. Circles 
and error bars are from Ref. 3. The labels of the horizontal axis 
represent k' = k + fco to reconcile the convention e^ = o = and 
the location of the minimum of the triplet dispersion in TlCuCl;? 
at momentum fc' = fco = (0, 0, 2-7r). 



0.006 




Fig. 2. Critical boson density as a function of temperature. 
Squares: magnetic field along the b direction. Tilted squares : 
magnetic field along the (1, 0, 2) direction (data from Oosawa et 
a/. 5 ). Full line: HFP result with the full dispersion relation. Dot- 
ted line: HFP result with the quadratic approximation for the 
dispersion relation ej. = fc 2 /(2m) and ks/m = 43.6 K. 



discrepancy at present. 

In the low-temperature limit, the quadratic approxi- 
mation would become asymptotically exact within the 
HFP theory, giving 8 ' 22 

1 /Tm\ 3/2 
n„(T ^ 0) = -C3/2 — (27) 



However, this <~ T 3 / 2 behavior (dotted line in Fig. 2) is 
only recovered at very low temperature and n cr (T) shows 
significant deviations from Eq. 27 already at 2 K. 

6. Critical field and interaction parameter U 

In the HFP approximation the critical field H C (T) is 
related to the critical density by 8 



(g/2) [H C (T) - H c (0)] = 2Un CI (T) 



(28) 
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Fig. 3. Critical field H c (normalized by the g factor and for two 
magnetic field direction: squares for H\\b and g = 2.06 and tilted 
squares for H J_ (1, 0, 2) and g = 2.23) as a function of the 
density n cr at the critical point (obtained from the the magne- 
tization m cr per dimcr by ra cr = m cr /(g/ia)). Data from Ref. 5. 
The full lines and the values of U and H c (0) are obtained from 
fits to Eq. 28. 



A linear relation between H C (T) and n cr (T) is indeed ob- 
served in the experimental data, as can be seen in Fig. 3. 
The least-square fits are performed in the low-density 
region (or cquivalcntly low-temperature). The values ob- 
tained for H c (0) are in good agreement with most of the 
previous estimates (see Tab. I). These fits also provide 
an estimate for U around 340 K. However, as it can bee 
seen in Fig. 4, a slightly smaller value for U (320 K) 
gives a critical field H C (T) which is in very good agree- 
ment with all the available experimental data, even at 
high temperatures. This value is close to that obtained 
from a similar HFP analysis (including a small magnetic 
exchange anisotropy) of the magnetization curves. 15 

In the literature the experimental data for H C (T) have 
been analyzed by fitting to the power-law (1). Values 
from (f) = 1.67 to 2.2 have been reported 5 ' 6 ' 8 ' 10 ' 11 and 
it has been suggested that the deviation from the HFP 
theory ((f) = 1.5) could be caused by fluctuations effects 
beyond the mean-field approximation. From our results 
it appears that a realistic dispersion relation 4 combined 
with an HFP treatment is able to reproduce the data ac- 
curately with a single adjustable parameter (U). It covers 
a wide temperature range from the very low temperature 
regime < 1 K where the quadratic approximation holds, 
up to - 8 K. 

7. Specific heat 

The specific heat of TICUCI3 under magnetic field was 
measured by Oosawa et al. w and shows a peak (with an 
asymmetric A shape) at the transition. In this section 
we compare these results with the prediction of the HFP 
theory. 

From Eq. 12 the expectation value of the energy per 
site in the normal phase (n c = 0) is 

(E) = -Un 2 + I -f±-i k f B (i k , T) (29) 



(2tt) 3 



The specific heat is obtained by differentiation with re- 



12 



10 



8 



6 



~i — 1 — 1 — 
Full e k U = 320K 
[g=2.06 and g=2.23] 
e k ~k 2 /(2m) [g=2.23] 

O Neutron H||b 

□ Mag. H||b 

A Heat H||b 

X Heat H||b 

X Heat Hl(l,0,2) 

O Mag. H±(l,0,2) 




Fig. 4. Critical field H C (T). Full line: HFP result with the full 
dispersion relation and U = 320 K (plotted for two values of the 
gyromagnctic factor). Dotted line: e k = k 2 /(2m) approximation 
(g = 2.23). Hexagons: Data from Ref. 6. Squares and crosses : 
Data from Ref. 10. Triangles and three-leg symbol: Data from 
Ref. 11. 



spect to temperature and we get: 



d 3 k 



de k 



(2tt)3 k 
On f d 3 k df B 



with ks — 1 and 
dn 



r d 3 k dfs 

J (2tt) 3 dT 



J (2ii)- i de k 

In the condensed phase, Eq. 12 gives 



(E) 



d 3 k 



f B (E k ,T) + - 



d 3 k 
(2n} 



(30) 



(31) 



-h + C 



After some algebra, we obtain the specific heat as: 

1 r d 3 k 
TJ J2n)- 



(32) 



dfB_\ 
dE ) 



+2U 



+ 



dn 
dT 

d 3 k e k 



n c — n — 



f d 3 l 

J (27T. 



YE k \^ + \ + Ed 4§ 



(33) 
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Fig. 5. Specific heat (per dimer) under an applied field (along 
the b axis) minus the specific heat in zero field. Full lines: HFP 
results with U = 320K. Circles: Measurements by Oosawa et 
al. 10 The results for the different values of H have been shifted 
by 0.04 for clarity. 

with 

l r d 3 k - dfB_ 
on _ T J (2^F efc dE 

df ~ i orr f tPfc £fc 7 a/ B M cff . !-\ 

(34) 

In Figs. 5 and 6 the HFP results above are compared 
with the data of Oosawa et al. for two magnetic field ori- 
entations. The theoretical curves reproduce qualitatively 
the A shape observed experimentally, although the height 
of the peak seems to be overestimated. 

8. Discussions 

In this paper, we have shown that taking the realis- 
tic dispersion relation determined from the microscopic 
theory and from the neutron scattering data, we can sig- 
nificantly improve the HFP approximation to explain the 
experimental data, especially the phase boundary curve 
H C (T). It is now evident that, in TICUCI3 the magnon 
dispersion curve is rather "steep" so that the quadratic 
approximation fails above a rather low temperature ~ 1 
K. 

It may be rather surprising that the HFP approxima- 
tion, which is generally believed to fail in the critical 
region, describes a wide range of experimental data pre- 
cisely. This appears to be the case, even though the HFP 
approximation still contains unsatisfactory features of 
predicting discontinuities in the magnetization and in the 



Fig. 6. Same as Fig. 5 with magnetic field H _L (1, 0, 2). 

specific heat at the transition. These discontinuities are 
considered to be an artifact of the HFP approximation. 
The true behavior of the magnetization in the model (2) 
is believed to be continuous and that of the specific heat 
to show a sharp cusp (negative exponent a, see Ref. 23) 
at the transition, which is classified as the 3-dimensional 
XY universality class. 

However, in fact, the experimental data on TlCuCl 3 
discussed in Section 7 does not show such a sharp sin- 
gularity and is rather similar to the HFP prediction. 
This may be explained by small anisotropies (breaking 
the U(l) symmetry around the magnetic field direction), 
which are expected to exist in any real magnetic sys- 
tem. The fact that the observed moment in the ordered 
(condensed) phase of TICUCI3 points to a constant direc- 
tion 6 suggest the presence of the anisotropy. Moreover, 
recently it is argued that a high-precision ESR measure- 
ment reveals the anisotropy. 24 Such anisotropies induce 
a small gap and should reduce the thermal fluctuations 
(and thus the specific heat) in the vicinity of the transi- 
tion, which could be also smeared out into a crossover. 
Since the breakdown of the HFP approximation is gen- 
erally due to the critical fluctuation, the reduction of 
the critical region caused by the magnetic anisotropies 
may actually make the agreement with the HFP predic- 
tions better, although we did not take any anisotropy 
into our calculation. Recently, an HFP calculation in- 
cluding a (small) magnetic anisotropy was carried out 
by Sirker et al. 15 and provided an improved description 
of the magnetization curves compared to that obtained 
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from the isotropic model. They also emphasized that the 
HFP approximation should be valid outside a narrow 
critical regime. 

Magnetic anisotropies are not the only corrections that 
may be added to the present model. Indeed, NMR mea- 
surements revealed that the transition to the ordered 
phase is (weakly) first order and accompanied by a si- 
multaneous lattice distortion 25 (see also Ref. 12). Spin- 
phonon interactions therefore seem to reduce the impor- 
tance critical fluctuations close to the transition while the 
resulting lattice distortion certainly induces some change 
in the magnetic exchange parameters. 26 An analysis of 
the consequences of such a magneto-elastic coupling is 
an interesting issue for further studies. 

Finally, let us comment on related theoretical works. 
Sherman et al. discussed that the agreement of the HFP 
result to the experiment is better if the "relativistic" 
form + A <~ Vc 2 fc 2 + A 2 is assumed for the magnon 
dispersion relation. 12 Our approach in this paper of mod- 
ifying the dispersion is actually the same to theirs. How- 
ever, we see no particular reason why we should take 
the relativistic form, although it may be a better ap- 
proximation for TICUCI3 than the quadratic one. In any 
case, ours would give a further improvement over Ref. 12 
within the HFP framework. 

In Rcfs. 13, 14, the phase boundary H C {T) is stud- 
ied numerically by a Monte Carlo method, for a dimcr 
system on a cubic lattice. The result should contain ef- 
fects from both the deviation of the dispersion from sim- 
ple quadratic, and the fluctuation beyond HFP While 
we cannot directly compare their result to ours as we 
deal with different models, the qualitative behavior is 
similar. Namely, they also observed the deviation from 
= 3/2 at higher temperatures, but the result seems to 
become closer to the = 3/2 as the temperature is low- 
ered. However, they suggest that this behavior including 
the deviation from = 3/2 could be universal and docs 
not depend on the particular dispersion, in a moderately 
low temperature regime. This is in contrast to our re- 
sult that the non-universal magnon dispersion explains 
the observed phase boundary H C (T) and its deviation 
from <\> = 3/2. The resolution is an open problem for 
the future. Numerical approaches would be also useful 
to clarify the effect of the (small) anisotropies. 
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Note Added 

After submission of this paper, Kawashima 27 clarified 
the question of the exponent <\> with numerical simu- 
lations of the 3D S=l/2 XXZ model as well as field- 
theoretical arguments. According to his results, in the 
limit of T — > 0, the HFP prediction = 3/2 is indeed ex- 
act. This is also consistent with our result that the phase 
boundary for a wide temperature range can be accounted 
within the HFP calculation using the realistic dispersion 
curve. 
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